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Abstract 

In order to clarify the influence of (the presence of) the broken time-reversal 
symmetry state (BTRSS) induced near the interface, tunneling conductance 
spectra in normal metal / d^2_y2-wave superconductor junctions are calcu- 
lated on the basis of the quasiclassical Green's function method. The spatial 
dependence of the pair potential in the superconductor side is determined 
self-consistently. We discuss two types of the symmetry on the BTRSS; i) 
(i^2_j^2-|-i,s-wave state and ii) d^2_y2-\-idxy-wave state. It is shown that the am- 
plitude of the subdominant component (is-wave or ida;j^-wave) is quite sensitive 
to the transmission coefficient of the junction. As the results, the splitting of 
the zero-bias conductance peak due to the BTRSS inducement is detectable 
only at junctions with small transmission coefficients for both cases. When 
the transmission coefficients are relatively large, the explicit peak splitting 
does not occur and the difference in the two cases appears in the height of 
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the zero-bias peaks. 

PACS numbers: 74.50.+r, 74.20.Rp, 74.72.-h 



Typeset using REVT^ 



I. INTRODUCTION 



Identifying of the pairing symmetry in high-T^ superconductors is important to clarify the 
mechanism of the superconductivity. A great deal of experimental and theoretical studies 
have revealed that the superconducting pair potential has (ia.2_y2-wave symmetry in the 
bulk stateB^. Since the pair potential of the (ij.2_j^2-wave superconductor is anisotropic, the 
amplitude of the pair potential near the surface or interface is significantly reduced. This 
suppression of the pair potential causes a very interesting phenomena, i.e. the formation of 
Andreev bound state (ABS) at the Fermi energy (zero-energy) near a specularly reflecting 
surfaceli when the angle between the lobe direction of the dx2-y2-wave pair potential and 
the normal to the interface is nonzero. This state is originated from the interference effect 
in the effective pair potential of the dx2-y2-wa:Ve symmetry through the reflection at the 
surface or interface. The ABS manifests itself as a sharp peak in the middle of the tunneling 
conductance spectra, the so-called zero-bias conductance peak (ZBCP),I and the consistency 
between theory and experiments has been checked in detailfflm. 

On the other hand, there still remains a controversial issue; formation of a broken 
time reversal symmetry state (BTRSS) at low temperature due to the mixing of subdom- 
inant or dxy-waye componentlisFiil as the imaginary part of the pair potential to 

the predominant dx^-y^-wave componentil. Theoretical studies based on the quasiclassical 
approximationii'ii and several lattice modelsilS reported the presence of the induced sub- 
dominant pair potential near the surface which breaks the time reversal symmetr yi. The 
resulting surface density of states with the BTRSS shows the splitting of the zero-energy 
peakil and the corresponding tunneling conductance shows the ZBCP splitting. It has also 
been clarified that the magnitude of the splitting depends on the induced subdominant pair 
potential near the surface. Actually, Covington et a/.H, Krupke et a/.ii, and Sharoni et a/.il 
reports the ZBCP splitting at low temperature and they ascribed the origin of the ZBCP 
splitting to the above BTRSS0. However, at the same time, there are many experiments 
which do not show the ZBCP splitting, and in these experiments, ZBCP survives even at 
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low temperaturesi'00llilli. Although there are several pre-existing theories which discuss 
the BTRSSH, almost all these theories treat the semi-infinite (i^2_„2-wave superconductor 
and the relevance to the actual experiments of tunneling spectroscopy has not been fully 
clarified yet. At this stage, it is an important problem to clarify the stability and the possible 
observability of the formation of BTRSS, i.e., the ZBCP splitting, for various condition of 
the junctions. There are several factors which determine the magnitude of the splitting of 
ZBCP; (i) inter-electron potential which induces the subdominant pair potential near the 
interface, (ii) transmission probability of the particles at the interface, (iii) orientation of 
the junctions, i.e., the angle between the normal to the interface and the crystal axis of 
(ia;2_y2-wave superconductor, (iv) finite temperature effect, (v) roughness of the interfaced, 
and (vi) impurity concentration in the superconductor^. 

Although there are several factors which determine the magnitude of the splitting of 
ZBCP, in this paper, we concentrate on (i), (ii), (iii), and (iv). First, we calculate the 
spatial variation of the pair potential in the normal metal / dx2_y2 -wave superconductor 
junction (N/D junction) in the presence of subdominant pair potential using a quasiclassical 
formalism for various conditions of the junctions. Using these results, we calculate tunneling 
conductance. 

The organization of this paper is as follows. In Sec. |T|, a formulation to calculate the 
spatial dependence of the pair potential, and the tunneling conductance is presented. In 
Sec. m, results of the numerical calculations for both (i^2_y2+is- and dx2_^y2+idxy-wave 



states are discussed in detail. Section IV is devoted to conclusions and future problems. 



II. THEORETICAL FORMULATION 

We study N/D junctions in the clean limit, where the normal metal is located at a; < 
and the (i^2_j^2-wave superconductor extends elsewhere. For the simplicity, two dimensional 
system is assumed and the x-axis is taken perpendicular to the fiat interface located at x = 0. 
When quasiparticles are in the xy-plaiae, a transmitted electron-like quasiparticle and hole- 
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like quasiparticle feel different effective pair potentials A(0_|.) and A(0_), with 0^. = and 
= 71 — (p. Here, (p is the azimuthal angle in the xy-plane given by [k^ + iky)/\k\ = e'*^. 
The barrier potential at the interface has a 5-functional form H6{x), where 6{x) and H are 
the 5-function and its amplitude, respectively. A cylindrical Fermi surface is assumed and 
the magnitude of the Fermi momentum and the effective mass are chosen to be equal both 
in the normal metal and in the superconductor. 

The quasiclassical Green's function method@'0i2h@ developed by Ashida et a/.ii'0 is 
used in order to determine the spatial variation of the pair potential self-consistently. In 
the following, we briefly summarize this quasiclassical method we used. We start with the 
Bogoliubov-de Gennes (BdG) equation for unconventional spin-singlet superconductorsHQ, 



EnUnir) = HoUnir) + J dr'A{r,r')v^{r'), 
EnVn{r) = -HoVnir) + J dr'A*{r,r')unir') 



2m 



fx, 



(1) 
(2) 

(3) 



where /i is the chemical potential, while Un{r) and Vn{r) denote the electron-like and hole-like 
components of the wave function. 



u„(r) ' 



V 



ifep-r 



(4) 



Vn{k,r) 

Here the quantities k and r stand for the unit vector of the wave number of the Cooper 
pair which is fixed on the Fermi surface {k = kp/\kp\), and the position of the center of 
mass of Cooper pair, respectively. After applying the quasiclassical approximation, the BdG 
equation is reduced to the Andreev equationiiiS, 



E„^'„(fc,r) 
A(k,r) 



ihvpk -V + A{k,r) f3\&„(fc,r) 



V 



A{k,r) 
-A*(k,r) 



(5) 
(6) 



where vp and fi{i = 1,2,3) stand for Fermi velocity and Pauli matrices, respectively. The 
wave function \[^„(fc,r) is obtained by neglecting the rapidly oscillating plane-wave part 
following the quasiclassical approximationilEl. The k dependence of A(fc, r) represents the 
symmetry of the pair potential. 

Now, we consider the case where a specularly reflecting surface or interface runs along 
the ^/-direction. In this case, the pair potential depends only on x since the system is 
homogeneous along the ?/-direction. It is convenient to introduce the following directional 
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Here ± represents the sign of the x component of the Fermi wave number kp^ and a{f3) = ±. 
We define a Green's function x, x') and a quasiclassical Green's function gai3{4>,x), 

$(°)(0,,x)$(f)t(0^,x) 



Gaf^{(f),x,x') = J2 



E„ 



ga/3{(p, X) ± i{%)ap = -2h\vF^\f3Gap{(l), X±0,x). 



(8) 
(9) 



In the above, 73 is the Pauli matrix in the directional spaceE^ and is the x-component 
of the Fermi velocity. The quasiclassical Green's function gai3{(f),x) obeys the Eilenberger 
equation,!! 



i\vFx\-^9ai3{(f),x) = - a 



iUmh + M4>a,x] 



9a(3[(P,X) 



+ 13 9afs{4>, X) \iuJmh + A{(f)fi, x) 

( \ 



A(0a,a;) 



V 



(10) 

(11) 



A(0„,x) 
-A*(0«,x) 

where uj^ is the Matsubara frequency. The quasiclassical Green's function can be written 
by the following evolution operator ?7a(0 



gal3{(l), X) = Ua{(f)a, X, x')gal3{4>, x)U ^ (0/3, X, x) 



(12) 



where Ua{4>a,x,x') satisfies the Andreev equation 



d 



-a 



(13) 



with Ua{(f)a,X,x) = 1. 



Considering a semi- infinite N/D junction geometry, the pair potential in the superconduc- 
tor side approaches to the bulk value A(0a, oo) at sufficiently large x. Hence, the evolution 
operator can be divided into a growing part and a decaying part: 



X, x') = Ai+)(0„, X, x')e«"(^-^') + Al-\<Pa, X, x')e-'''^(^-^'), 
1 



-<i>i+\(Po.,xf<^i-\4>^,x')f2, 



, y<ra,X,x') = — $i \(l)a,xf^l^\(l)^,x')f2, 



(14) 



(15) 



where 



\VFi 



|A(0„, oo) 



= T<l>W(0„,x)f2<l'i-)(0.,x) = -T<-)(0,,x)f2<l>W(0«,a;) 
= constant. 



In the above, '^$^"^(0q,, x) denotes the transposition of $^"^(0a,a;) 



(16) 



Retaining the most divergent term in semi-infinite limit, we find the quasiclassical Green's 
function gaai4>yx) in the superconductor side given bj''^^''^ 

2Asa{x) 
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.X = 1 



M^Saix)] 



(17) 



where 



As+{x) = Ai \<j)+, X, L)aL+^ L, O)^Nf/+(0+, 0, x)e 



As+(x, 0)RTS(U+{(f>+, 0, x) 



As+(a;,0) oc 



^ ui H(j)4^, x) ^ 



4 ^(0+,a;) 



(18) 
(19) 



In the above, the matrix represents resistance at the interface which is given byil 



1 

1 - M<P) 



crN(0) 



4 cos 0^ 

4COS02 + Z2' 



(20) 



where o"n(0) stands for tunneling conductance when the system is in the normal statee^ and 
Z is the effective barrier height at the interface with Z = 2mH / {fi^k-p)- In order to obtain 
the quantity ?7+(0+,O,x) in eq.(^8|), we rewrite y4g+(x) asS 



(21) 



where Va{x) = {—i)v!^\(f)a, x) /u^^^cpa, x) and J-'+{x) = iX+{x) /Y+{x) obey the following 
Riccati type equations 
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^|fFx|7^I'a(a;) =a -2uJm'Do,{x) + A{(f)a, x)Vl{x) - A*(0q, 
ax L 

d 

Mvf^\—J^+{x) = 2uJrnT^{x) + iS!" {(^^ , x)T\{x) - A(0+,a;). 
ax 



X 



(22) 
(23) 



We can write the quasiclassical Green's function in a compact formil, 

\ 

il?+(x) -I?+(x)^+(x) 



,x) = i 



{ 1 iJ^+(x) 



l-I?+(x)n(x) 



V 



(24) 



Initial conditions of these equations are 

VJoo) 



Urn + Oi^a 

The pair potential is given byHSSiliiS 



n(o) 



1 - O-n(0) 

D_(o) • 



(25) 



1 ^ ^ 



(26) 



with ^ {(j)^,x) = — ^++(— 0+,x), where c^c is the cutoff energy and [gaa{(pa,x)]i2 means 

the 12 element of gaa{4'a,x). Here V{(j),(f)a) is the effective inter-electron potential of the 
Cooper pair. In our numerical calculations, new A(0q,,x) and gaa{(pa,x) are obtained using 



eqs.(^)-(|2J) and eq. (^61) . We repeat this iteration process until the sufficient convergence 
is obtained. 
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Next, we calculate the tunneling conductance spectra based on the self-consistently de- 
termined pair potential. The resulting normalized tunneling conductance a^ieV) with the 

bias voltage V is given byillii 

r-VS , foo ^ /E -eV 



/n/2 i-oo / — \ 

dcj) / dEa^{(l))as{E, 0)sech^ cos< 
, w , -■n'/2 J -co \ Zk-qI / /ni-r\ 

^T{eV) = '-—j-^ — fE-eV\ ' 

/ d(j) rfEaN(0)sech^ cos< 

i + aN(0)|rs+(i^,0+,o)P + kN(0)-i]|rs+(i^,0+,o)nrs4i^,0-,o)P ..8^ 
""'^ ''^^ |i + K(0)-i]rs+(i?,0+,o)rs_(i^,0-,o)|2 ' ^ ^ 

In the above, rsa(-E, 0a, a;) is obtained by solving following equations, 

\ri\v^^\^Y^^{E, 0+, X) = 2ETs+iE, 0+, x) - A(0+, x)Tl4E, 0+, x) - A*(0+, x), (29) 
ih\vF,\^Ts-{E, 0-, x) = 2Ers_(S, 0-, x) - A*(0_, a;)r2_(E, x) - A(0_, x). (30) 

In the following, almost calculations are performed on the temperature T/Tc = 0.05, where 
Tc is the critical temperature of the bulk (i^.2_y2-wave superconductor. 



III. BROKEN TIME REVERSAL SYMMETRY STATE NEAR AN INTERFACE 
OF A L>^2 y2-WAVE SUPERCONDUCTOR 

In this section, the spatial dependence of the self-consistently determined pair potential 
and the corresponding tunneling conductance are presented for the dx2_y2-wave supercon- 
ducting state. In the middle of the c?a,2_j^2-wave superconductor, the pair potential is given 
by A(</)q, oo) = Aq cos[2(</) — 6)], where 6 is the angle between normal to the interface and 
the lobe direction of the dx^^yi-wave pair potential, i.e., the angle between the x-axis and 
the crystal a-axis of the dx2-y2-'wave. In this paper, we choose various ^ (0 < < 7r/4) by 
changing the magnitude of Z and (T^^^). 



A. dx2_y2+ is- wave state 



In this subsection, we show the spatial dependence of the pair potential and the resulting 
tunneling conductance of the (ij.2_j^2+is-wave state realized near the interface of the N/D 
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junction. The spatial dependence of the pair potential is expressed as 

A(0, x) = A,(x) cos[2((/) - e)] + A,(a;), (31) 

where Ad(x) and As(a:) correspond to the amplitude of the (ia;2_y2-wave and s-wave super- 
conducting states, respectively. The attractive potential V^(0, (p') is given by 

V{<P, 0') = 2Vd cos[2(0 - e)] cos[2(0' - 9)] + (32) 

where Vd and K denote the attractive potential of predominant (ia;2_y2-wave and subdomi- 
nant s-wave, respectively, and they are given as 

Vd = J , (33) 

= ^ ^. (34) 

log — + 711^ 

Ts 1^ rpTn + 1/2 

Here, Ts denotes the transition temperature of s-wave component of the pair potential 
without predominant (ia,2_y2-wave component. For 9 = {] oi 9 = 7r/4, only Re[Arf(x)] 
and Im[A5(x)] are nonzero. The spatial dependence of the pair potentials Re[Ad(x)] and 
Im[As(a;)] is plotted in Fig. |l|(a) for various 9 with T^jTd = 0.2 and Z = 3. The x-axis of 
Fig. |l|(a) is normalized by = ^"^f / Aq which is the coherence length of the superconductor. 
For ^ = [(100) surface], since Im[As(x)] = is satisfied, the time reversal symmetry is not 
broken and the amplitude of Re[Arf(x)] is not suppressed at the interface. By changing the 
angle 9 from zero, the magnitude of Re[A(i(a;)] is reduced near the interface, while Im[As(x)] 
is induced at the interface@. The suppression of Re[A(i(a;)] is originated from a depairing 
effect that the effective pair potentials A(0_|_,O) and A(0„,O) have reversed contribution to 
the pairing interaction for certain range of for 9^0. When Re[Ad(x)] is suppressed at 
the interface, the quasiparticle forms the ABS near the interface at zero-energ yi. The ABS 
is unstable with the introduction of s-wave attractive potential, then Im[A5(x)] is induced 
at the interfaceii. The magnitude of Im[As(x)] becomes maximum at ^ = 7r/4, where the 
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above suppression effect is most significant. In Fig. |l](b), Re[Ad(a;)] and Im[As(x)] are plot- 
ted for various Z with T^jT^ = 0.2 and 6 = n/A. Even if tlie BTRSS becomes to be most 
stable at 6 = 7r/4, when the height of barrier is small, the magnitude of the subdominant 
imaginary component of As{x) is not induced at all. The induced imaginary component of 
A<j(x) is enhanced with the increase of Z. 

The spatial dependence of the pair potentials near the interface with the intermediate 
angle {6 = vr/6) is shown in Fig. ^ for various height of barrier. In such a case, both 
Im[Arf(x)] and Re[As(x)] becomes nonzero and the spatial dependence is much more complex 
as compared to that for 6 = ot 6 = n/A. The amplitudes of Im[Ad(a:;)], Re[As(x)], and 
Im[As(x)] are enhanced for larger magnitude of Z, where the suppression of Re[A(;(x)] is 
significant. However, the amplitudes of Im[Arf(a;)] and Re[As(x)] are one order smaller than 
that of Im[As(x)]. 

Next, we look at the magnitude of subdominant components of the pair potential at 
the interface, Im[As(0)], Im[A(i(0)], and Re[As(0)], for various T^, Z, and 6. As shown in 
Fig. ^(a), the magnitude of Im[As(0)] increases monotonically with for fixed 6 and Z, 
and it is enhanced for larger magnitude of Z. In other words, the amplitude of Im[As(0)] is 
sensitive to the transmission probability of the junctions. In Fig. ^(b), Im[As(0)] is plotted 
as a function of 6 for sufficiently larger magnitude of Z{= 5.0). For ^ = 0, i.e., junction 
with (100) interface, the magnitude of Im[As(0)] is negligibly small near the interface even 
at the larger magnitude of T^. The magnitude of Im[As(0)] is a monotonically increasing 
function with the increase of 6 and has a maximum at = 7r/4. As seen from Figs, ^(c) and 
^(d), both the magnitude of Im[Ad(0)] and Re[As(0)] is enhanced and has a maximum at a 
certain 6. In the intermediate 6, i.e., 6 ^ or 6 ^ 7r/4, the magnitude of A(0_(.,x) and that 
of A (</)„, x) does not coincide any more, the interference with the quasiparticle and the pair 
potential becomes complex. Then, not only Im[A<,(0)] but also Re[As(0)] and Im[Ad(0)] 
become nonzero. 

Using self-consistently determined pair potentials, let us look at the normalized tunneling 
conductance aT^eV). In order to clarify the temperature T dependence of aT^eV), we choose 
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T = in the left panels of Fig. ^ and T = O.OSTc in the right panels. Only for 6 = 0, line 
shape of aT^eV) is similar to that of the bulk density of states of (i^2_y2-wave superconducting 
state. In other cases, aT{eV) has a zero bias enhanced line shape. As clarified in previous 
literatures^, when 6 deviates from zero, since the injected and reflected quasiparticles have 
a chance to feel the sign change of the pair potentials, zero-energy ABS is formed at the 
interface. This zero-energy ABS causes the ZBCP when the magnitude of Tg is small. With 
the increase of the magnitude of T<j, the zero energy ABS is unstable and s-wave subdominant 
component is induced which breaks time reversal symmetry and it blocks the motion of the 
quasiparticles. Then, the energy levels of bound state shift from zero and the local density 
of states has a zero-energy peak splitting. The resulting aT^eV) has a ZBCP sphtting as 
shown in Fig. §(b). However, with the increase of T, the slight splitting of ZBCP fades out 
due to smearing effect by finite temperature and the resulting aT^eV) has a rather broad 
ZBCP [see Fig. |(c) and Fig. |(d)]. 

Next, we concentrate on how a^ieV) is influenced by the transmission probability of the 
junctions, i.e., the magnitude of Z. In Fig. |^, aT^eV) with 9 = 7r/4 is plotted for T = (left 
panels) and T = O.OSTc (right panels). For the junctions with high transmissivity, a^^leV) 
has a ZBCP, [see Fig. ^a)] and the magnitude of o"t(0) is firstly enhanced with the increase 
of Z. In this case, the predominant drc2_y2-wa:Ve component only exists near the interface as 
shown in Fig. |l](b). However, with the increase of Z, a^leV) starts to have a ZBCP splitting 
at a certain value of Z, where the magnitude of subdominant component Im[As(x)] at the 
interface becomes the same order as that of the predominant component Re[A(;(x)]. For 
sufficiently larger magnitude of Z, aT^eV) has a ZBCP splitting, [see Fig. ^(b) and |(c),] 
and the magnitude of o"t(0) decreases with the increase of Z. However, the above obtained 
results are influenced by flnite temperature effect. The right panels of Fig. ^ is shown 
for the tunneling conductance in T/T^ = 0.05. The slight enhanced structure of aT{eV) 
at eV = ±Ao in Fig. |^(a), ||(b), and |^(c) is invisible due to the smearing effect by flnite 
temperature [see Fig. ^c), |§(d), and|^(e)]. With the increase of the magnitude of Z, aT^eV) 
has a ZBCP with tiny dip even at Z = 2.5, where the order of the amplitude of Im[As(0)] is 
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0.2Ao. With the further increase of aT^eV) has a ZBCP sphtting [see Fig. ^(f)], however 
the degree of the sphtting is significantly weakened as compared to the corresponding curves 
in Fig. |(c). 

FinaUy, we look at the relation between the position of the splitted peak and the mag- 
nitude of Im[As(0)]. In Fig. ^(a), aT^eV) is plotted for various Tg with Z = 5 and T = 0. 
As shown in Fig. |](b), the magnitude of the induced subdominant imaginary component of 
Im[A,(0)] is about O.I6A0, O.SAq, and 0.42Ao for T./T^ = 0.1, T./T^ = 0.3 and T./Td = 0.5, 
respectively. The corresponding aT^eV) has a splitted peak locating at ±0.16Ao, ±0.3Ao, 
and ±0.42Ao, respectively. With the increase of T, the height of these peaks are drastically 
suppressed as shown in Fig. ^b). Finally, we show how the line shape of a^ieV) changes 
at the temperature T = Tg {Tg = 0.12Tc) where Im[As(x)] becomes nonzero. As seen from 
Fig. H(c), the magnitude of crr(O) is reduced with the introduction of Im[As(x)]. 

At the end of this subsection, we can summarize that even in the presence of BTRSS, the 
resulting aT{eV) does not always have a clear ZBCP splitting due to the finite temperature 
effect when the magnitude of the transmission probability of the junctions is not low. 

wave state 

In this subsection, we study spatial dependence of the pair potentials of the dx2_y2+idxy- 
wave state and the resulting tunneling conductance in N/D junctions. The pair potential is 
given by 

A(0,x) = A,(a;)cos[2(0-^)] + Arf^^(x)sin[2(0-^)], (35) 

where A^^^^ (x) is an amplitude of the rf^jy-wave superconducting state and a complex number. 
The attractive inter-electron potential V{(f), (f)') is given by 

V{(t>, <P') = IVd cos[2(0 - e)\ cos[2(0' - e)\ + 2V^^^ sin[2(0 - e)\ sin[2(0' - 0)], (36) 

where Vd and Vd^y stand for the attractive potential of predominant (i^2_j^2-wave and sub- 
dominant (i^y-wave, respectively, and they are given as 
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Vd = '-^ ^, (37) 

2Tlk-aT , ^ 

Vdxy = ^ 1 . (38) 

The spatial dependence of the pair potentials with a finite transmissivity for Z = 3.0 and 
^ = 7r/4 is shown in Fig. 0(a). As in the case for (i^2_j^2+is-wave state, the amplitude of 
Im[Afj^j^ (x)] vanishes at 6' = 0. For 6 = tt/A, the suppression of the magnitude of Re[A(i(a;)] is 
most significant, while Ad^y{x) is induced at the interface. At this 6, the d^y-wave component 
is not affected by depairing effect seriously since the lobe direction of dxy-v/ave pair potential 
is parallel or perpendicular to the interface as in the case for predominant d^2_y2-w8ive with 
6 = 0. As shown in Fig. 0(b), the amplitude of Im[Ad^j^(x)] is enhanced with the increase of 
Z. 

The spatial dependencies of Re [A^^ (a;)], Im[Arf(x)], Re[As(x)], and Im[As(a:;)] are plotted 
in Figs. H(a), ^(h), ^{c), and||(d), respectively, with the intermediate 9 {9 = n/6) for various 
Z. In such a case, the spatial dependence is much more complex as compared to that for 
6' = or = 7r/4, and both Im[Ad(a;)] and Re[Ad^j^(x)] become nonzero. The amplitudes 
of Im[Arf(x)], Re[Arf^j^(x)], and Im[As(x)] are enhanced for larger magnitude of Z, where 
the suppression of Re[Arf(a;)] is significant. The remarkable feature is that the amplitude of 
Im[Ad(a;)] and Re[Ad^y(a:;)] can become the same order of Im[As(a;)] for larger Z. 

Next, we look at the magnitude of subdominant components of the pair potential at 
the interface, Im[Ad^y(0)], Im[Ad(0)], and Re[Ad^j^(0)] for various parameters Td^y, Z and 9 
shown in Fig. ^ The magnitude of Im[Arf^j^(0)] increases monotonically with Td^y for fixed 
9 and Z. It is also enhanced for larger magnitude of Z. Comparing the corresponding 
situation of dx2_y2 + is- wave shown in Fig. ^(a), the magnitude of Im[Ad^^(0)] is suppressed 
for the small magnitude of Td^y. In Fig. ^(b), Im[Ad^j^(0)] is plotted as a function of 9 for 
sufficiently larger magnitude of Z{= 5.0). The magnitude of Im[Ad^j^(0)] is a monotonically 
increasing function with the increase of 9 and has a maximum at = 7r/4. The amplitudes of 
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Im[Ad(0)] and Re[A(i^j^(0)] are negligibly small for 6 = and 6 = it /A, and has a maximum 
at a certain 9 as shown in Figs, ^(c) and ^(d), respectively. 

The resulting normalized tunneling conductance aT^eV) is plotted in Figs, ^{a) and 
|10|(b) with T = and the corresponding quantities with T = O.OSTc is plotted in Figs. p]0|(c) 
and |Tn|(d) for various 6. Only for = 0, line shape of a^^leV) is similar to that of the bulk 
density of states of dx2_y2-wa:Ve superconducting state. In other cases, aT{eV) has a zero-bias 
enhanced line shape due to the formation of ABS [see Fig. §(a)]. The ABS causes the ZBCP 
when the magnitude of T^^^ is small. However, with the increase of the magnitude of T^^^, the 
zero-energy ABS is unstable and subdominant ci^^y-wave component is induced. The energy 
levels of bound state shifts from zero and the local density of states has a zero-energy peak 
splitting. The resulting aT{eV) has a ZBCP splitting as shown in Fig. |10](b). As compared 
to the case of dx2^y2+is state, the magnitude of o"t(0) is not reduced seriously since the 
effective dxy-wave pair potential felt by quasiparticles is distributed from — Im[Ad^j^(0)] to 
Im[Arf^j^ (0)]. With the increase of temperature, the ZBCP splitting is not visible any more, 
[see Fig. |rU|(c) and Fig. |ig(d)]. 

Next, we concentrate on how a^leV) is influenced by the transmissivity of the junctions, 
i.e., the magnitude of Z. In Fig. |n|, a^leV) with 6 = 7r/4 is plotted for T = (left 
panels) and T = O.OSTc (right panels). As shown in Fig. |ll], aT{eV) with 6 = n/A depends 
on the transmissivity of the junction crucially. For the junctions with high transmissivity, 
aT^eV) has a ZBCP [see Fig. |ll](a)] and the magnitude of <Jt{0) is firstly enhanced with 
the increase of Z. In this case, only the predominant (i^2_j^2-wave component exists near 
the interface as shown in Fig. |l](b). However, with the increase of Z, a^leV) starts to have 
a peak splitting at a certain value of Z, where the magnitude of the induced subdominant 
component at the interface Im[A(i^j^(0)] overlaps that of the predominant one Re[A£;(a;)]. For 
sufficiently larger magnitude of Z, anci^eV) has a ZBCP splitting [see Figs. |ll](b) and |ll](c)] 
and the magnitude of crT(O) decreases with the increase of Z. The right panels of Fig. ^ 
are shown for the a^leV) with finite temperature. Only for the junctions with the lowest 
transmissivity {Z = 5), the subtle splitting of ZBCP appears in aT^eV). For T = O.OSTc, 
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similar line shapes of arieV) with smeared ZBCP splitting are obtained [see Fig. 0(d), 
Fig. |ll](e) and Fig. |ll](f)]. However, the slight enhanced structure of axieV) at eV = ±Aq 
vanishes. 

Finally, we look at the relation between the position of the splitted peak and the magni- 
tude of Im[A(i^^(0)]. In Fig. |T^(a), we show T^^^ dependence of the tunneling conductance 
spectra aT^eV) for 6* = with Z = 5. As shown in Fig. |(b), the magnitude of the induced 
subdominant imaginary component of Im[A(i^j^(0)] is about O.ISAq, 0.35Ao, and 0.52Ao for 
Td,y/Td = 0.1, Td.jTd = 0.3 and Ta^jTd = 0.5, respectively. The corresponding aT{eV) 
has a splitted peak locating at ±0.16Ao, ±0.3Ao, and ±0.42Ao, respectively. The width of 
the two splitted peaks depends on the temperature as shown in Fig. |l3(b). Next, we show 
how the line shape of aT^eV) changes at the temperature T = T^-^ [Td~y = 0.07Tc] where 
Im[Ad^j^(x)] becomes nonzero. As seen from Fig. |l3(c), the magnitude of o"t(0) is reduced 
with the introduction of Im[As(x)]. 

At the end of this section, we can summarize that even in the presence of BTRSS, the 
resulting a^leV) does not always have a clear ZBCP splitting due to the finite temperature 
effect when the transmission probability of the particles at the interface is not low. In 
order to observe the ZBCP splitting which is one of the evidence supporting BTRSS, we 
must measure a^leV) for the junctions with low transmissivity with (110) oriented interface 
{6 = 7r/4) at low temperatures. In such a case, we can classify the the d^2_y2+is- and 
dx2_y2+idxy-wave state through the magnitude of o"t(0) and the width of the ZBCP splitting 
for the junctions with low transmission probability. 



IV. CONCLUSION 

In this paper, spatial dependence of the pair potential in the N/D junctions is determined 
on the basis of the quasiclassical theory in the presence of subdominant component of the 
pair potential near the interface. We clarified the influence of the spatial variation of the 
pair potentials on the tunneling conductance spectra for various conditions of the junctions. 
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We selected two kinds of subdominant components s- and dxy-wave which are induced as 
dx2_y2+is- and (i^2„j^2+i(ia,y-wave state, respectively. The amplitude of As{x) [Arf^^(x)] is 
enhanced with the increase of the magnitude of Z, Ts[Td^y], and (0 < < 7r/4). In the 
intermediate 9, both Im[Ad(x)] and Re[As(a;)] {Re[Ad^,^]} becomes nonzero, although these 
magnitudes are small as compared to those of Re[A(i(x)] and Im[As(x)] {Im[Arf^j^]}. The 
resulting aT^eV) depends on Z, Tg (or T^^y), and 0. For fixed Z, the magnitude of o"t(0) 
increases with the increase of (0 < < 7r/4) for small magnitude of Ts/Td \Td^y/Td\. 
While crT(O) increases and decreases again with the increase of ^ (0 < < 7r/4) due to the 
formation of d^^-yi+is- [c?^2_y2+ic?2:j,]-wave state for sufficiently larger magnitude of Ts/Td 
[Td^y/Td]. For fixed 9 and T^/Td [Td^y/Td], anc{eV) has a ZBCP for small magnitude of Z, 
while it has a ZBCP splitting for sufficiently larger magnitude of Z. Using junctions with 
small transmissivity, we can distinguish (ij.2_y2+is-wave state from dx2_y2-\-idxy-wa:Ve state 
since the magnitude of o"t(0) for 6/^:2 _y2+is-wave state is much more reduced as compared 
to that for (i^.2_y2+i(i^j^-wave state as seen from Fig. ^ [Fig. ^ to Fig. |10| [Fig. |12[. By taking 
into account finite temperature effect, the degree of the ZBCP splitting is suppressed and 
the fine structure at eV ~ Aq in aT^eV) becomes invisible. 

In the light of our theory, one of the reason for the absence of ZBCP splitting for many 
tunneling experiments may be due to their high transmissivity of the junctions (small mag- 
nitude of Z) and high temperatures. If we choose as 90K, unless the transmissivity of the 
junction is sufficiently low, it is difficult to observe ZBCP splitting clearly at 4.5K. In order 
to see the ZBCP splitting clearly, we must observe much more low temperature (T < 4.5K) 
using a clean junctions with low transmissivity for 9 = it /A. 

In this paper, we have chosen a free electron model with cylindrical Fermi surfaces. In 
order to compare actual tunneling conductance, we must calculate spatial dependence of the 
pair potential and the tunneling conductance by taking into account of the actual shape of 
Fermi surface. For this purpose, it is a promising way to perform the calculation based on the 
t-J model with Gutzwiller approximation@, where the doping dependence is naturally taken 
into account. Using this model, we succeeded to explain detailed line shape of the scanning 
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tunneling conductance spectra around Zn impurityii^0. As regards quasiparticle states : 
the interface, the pre-existing calculations are performed for infinite barrier limitS'i'i. In 
order to compare actual tunneling experiments, we must calculate the N/D junctions based 
on the t-J model using Gutzwiller approximation. In the t-J model, the subdominant dxy- 
wave component is hard to be realized^ since a nearest neighbor attractive potential is taken 
into account. In order to discuss the (i^^y-wave component, we must consider much more long 
range interaction including three site hopping term. 

It is also revealed by recent theory that the BTRSS influences significantly on the Joseph- 
son effect&S. The previous theory assumes a sufficient larger magnitude of barrier^. In 
order to compare the existing experiments, we must calculate Josephson current in the 
presence of BTRSS for arbitrary transmissivity of the junctions. 
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FIGURES 




X/ ^0 ^^^0 



FIG. 1. Spatial dependences of the pair potential at Tg/Td = 0.2; (a) near the interface for the 
angle 6 between x-axis and crystal axis at Z = 3.0, and (b) near the (110) interface [6 = 7r/4] for 
various height of barrier Z. T = O.OSTc. 
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FIG. 2. Spatial dependence of the pair potentials near the interface with = 7r/6 and 
Ts/T^ = 0.3; (a) real part of Arf(a;), (b) imaginary part of A(i(a;), (c) real part of As(a;), and 
(d) imaginary part of As(x). T = 0.02Tc. 
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FIG. 4. The normalized tunneling conductance in the dj.2_j^2+is-wave state near the interface 
at Z = 3.0. Ts/Td = 0.05 for (a) and (c). Ts/Ta = 0.1 for (b) and (d). T = for (a) and (b). 
T/Tc = 0.05 for (c) and (d). 
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FIG. 5. Tunneling conductance for dj.2_j^2+is-wave state with TgjT^ = 0.2 and 6 = 7r/4 for 
various Z; (a) [(d)] low barrier {Z = 1.0, 2.0) (b) [(e)] middle barrier {Z = 2.5), and (c) [(f)] high 
barier (Z = 3.0, 5.0). T = for (a), (b), and (c). T/Tc = 0.05 for (d), (e), and (f). 
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FIG. 7. Spatial dependences of the pair potential with Td^^/Ta = 0.2. T = O.OSTc. (a) Z = 3.0 
for various 9. (b) ^ = for various Z. 
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FIG. 8. Spatial dependence of the pair potentials near the interface with 6 = tt/G and 
Td^y/Td = 0.3; (a) real part of Ad{x), (b) imaginary part of Ad{x), (c) real part of Ad^y{x), 
and (d) imaginary part of Ad {x). T = 0.02rc. 
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FIG. 10. The normalized tunneling conductance for d^2_y2+idxy-waYe state with Z = 3.0. (a) 
Td,y/Td = 0.1 and T = 0.0. (b) Ta^jTa = 0.2 and T = 0.0. (c) T^^jTa = 0.1 and T = O.OSTe. (d) 
Td../rd = 0.2andT = 0.05r,. 
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FIG. 11. Tunneling conductance for dj.2_j^2+i(i-rj,-wave state with = tt/A and T^^^/Td = 0.2 
for various Z. (a) [(d)] low barrier {Z = 1.0, 2.0) (b) [(e)] middle barrier {Z = 2.5), and (c) [(f)] 
high barier {Z = 3.0, 5.0). T = for (a), (b), and (c). T = O.OST^ for (d), (e), and (f). 
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